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ABSTRACT 

Given the number of recently discovered galaxy-galaxy lens systems, we anticipate that a gravita- 
tionally lensed supernova will be observed within the next few years. We explore the possibility that 
stars in the lens galaxy will produce observable microlensing fluctuations in lensed supernova light 
curves. For typical parameters, we predict that ^70% of lensed SNe will show microlensing fluctua- 
tions > 0.5 mag, while ^25% will have fluctuations > 1 mag. Thus microlensing of lensed supernova 
will be both ubiquitous and observable. Additionally, we show that microlensing fluctuations will 
complicate measurements of time delays from multiply imaged supernovae: time delays accurate to 
better than a few days will be difficult to obtain. We also consider prospects for extracting the lens 
galaxy's stellar mass fraction and mass function from microlensing fluctuations via a new statistical 
measure, the time-weighted light curve derivative. 

Subject headings: gravitational lensing - galaxies: stellar content - supernovae: general 
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1. INTRODUCTION 

Within the past year there has been an explosion 
in the number of known galaxy-galaxy strong lens 
systems due to t he success of the Sloan Lens ACS 
(SLAGS') survey lIRolton et aU 120041 WM. iTreu et a,U 
120061 iKoopmans et al.ll2006D . The lens galaxies consist 
of a subset of the Sloan Digital Sk y Survey (SPSS) Lumi- 
nous Red Galaxy (LRG) catalog l)Eisenstein et al.ll2001|) 
and were initially selected based on the presence of sev- 
eral "rogue" oxygen emissi on lines in the othe rwise well 
characterized LRG spectra ijBolton et al.ll200"4j) . The hy- 
pothesis that these emission features originated from an 
unresolved background galaxy within the 3" SDSS fiber 
suggested a substantial probability of lensing. Follow-up 
HST snapshots with the ACS camera reveal ed that many 
of the sour ce galaxies are in fact lensed l)Bolton et al.l 
l200H 12001 . 

Oxygen emission lines tra ce star formation in both the 
local and distant universe l|Kennicuttlll992t iHogg et all 



f998), so it seems likely that the SLACS sources are un- 
dergoing intense star formation. Coupled with the sheer 
number of lenses discovered with SLACS (~30 at present 
count), this opens up the exciting possibility of observing 
lensed supernovae (SNe) in the near future. 

Lensed SNe offer a variety of inter e sting applica- 
tions. For example, iQguri fc Kawanol l)2003l) showed 
that SN magnification distributions can be used to break 
the degeneracy between the density profile of the lens 
galaxy and the me asured value of the Hubble constant. 
iCoorav et all (|2006) investigated the possibility of using 
weakly lensed SNe to ma p out large-scale density fluc- 
tuations in the universe. iKolatt fc Bartelmannl (1998) 
proposed the use of SN magnification distributions to 
eliminate the mass-sheet degeneracy in cluster mass re- 
construction. In this paper, we point out that strongly 
lensed SNe will also be microlensed by the stars in the 
lens galaxy, and demonstrate that the microlensing fluc- 
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tuations will contain valuable information about the stel- 
lar population of the lens galaxy. 

The first lensed SN to be discovered will probably be 
Type II, as opposed to the Type la SN e used to measure 
the expansion history of the universe (lRiess_et "alli"998t 
Perlmutter et all999tlTonrv et al2003tlKnop et al2003 : 



Riess et alll2004j) . While SNe II appear to be somewhat 
less common than SNe la when averaged over galaxies of 
all Hubble types, SNe I I have a higher rate among spiral 
and i rregular galaxies (JKennicuttJ 119921: IMannucci et al 
HHJ. The presence of [O II] emission lines in the SLACS 
sources suggests that they are predominantly gas-rich 
spirals and irregulars. Since star formi ng galaxies typ- 
ically host one SN II every 50-100 years ijMannucci et al 
2005) , we expect a SN II to be observed in one of the ~30 
SLACS source galaxies within the next several years. Al- 
though there is scatter among the peak brightnesses of 
SN IPs, the mean B-band luminosity is approximately 
M B ~ -17.5 llYoung fc Branch! ll989|K This implies an 
apparent magnitude of a SN in the SLACS source popula- 
tion (which has a mean redshift of z s « 0.53) of m ~ 24.5. 
Since there is an additional amplification /io ~ 10 due to 
the main lens galaxy, the apparent magnitude for multi- 
ply imaged SNe in the SLACS sample should be m ~ 22, 
a value that is certainly within reach of existing tele- 
scopes with reasonable exposure times. 

In the SLACS sample, discoveries of lensed SNe will 
come from monitoring known lens systems. The com- 
plementary approach is to mine a catalog of observed 
SNe to find those that are lensed. This approach could 
be used with the Sup ernova Acceleration Probe (SNAP; 
lAldering et al.l 12004 also http://snap.lbl.gov), which is 
intended to observe ^2000-3000 SNe la per y ear down to 
26-27th magnitude. IQguri. Suto. fc Turner! ((2003) fore- 
cast that SNAP will discover 2-20 lensed SNe la per year, 
with the uncertainties due mainly to the star formation 
rate at high redshift. Although SNAP will focus on SNe 
la, the deep imaging will naturally reveal many SNe II as 
well, raising the total number of lensed SNe by a factor 
of -2. 

While there are many categories a nd subcategories of 
observed SNe (see iFilippenka 119971 for a thorough re- 
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view) , the physics behind SN explosions has been studied 
extensively and both light curves and s pectra are rea- 
sonably well reproduced in simulations l | Eastman et al 
119941: lEastman. Schmidt, fc Kir shnerl 11996(1 . The pro- 
genitors of SNe II ar e thought to be massive (> 8-10 M Q ; 
iLeonard et all 120011) stars that undergo core collapse 
when the amount of energy released in fusion reactions is 
no longer sufficient to counteract self-gravitational forces. 
The resulting explosion creates an expanding envelope 
that can be modeled with the E xpanding Photosphere 
Method (Kirs hner fc Kwarall974(l . Light curves of some 
SN II have been observed with time sampling a pproach- 
ing 1 day (e.g., SN1999em: [Leonard et a1l2001l) . 

Light that reaches us from a distant SN will of course 
be affected by all intervening matter, such as dust and 
gas, leading to photometric and spectroscopic differences 
between the intrinsic and observed light. Lensed SNe are 
unique in that they may also exhibit gravitationally in- 
duced variability due to intervening stars lying within 
the lens galaxy. This phenomenon known as microlens- 
ing is familiar from quasar l ensing. with the best known 
example being Q2237+030 ijWozniak et all 120001) . The 
brightnesses of the lensed images fluctuate on time scales 
of months to years due to relative motion of the lens 
galaxy and source quasar. For a d etailed example o f 
microlensing light curve analysis see iKochanekl (J2004), 
and for detailed treatments of quasar microlensing theory 
see iWvithe fc Turner! (|2001fl: ISchechter fc Wambsganssl 
( 20021: ISchechter et a,U l|200i) . 

In a microlensed SN, by contrast, what changes with 
time is the size of the source. 3 The eff ects of source size 
in mi c rolensing have been studied (e.g., [ Mortonson et alJ 
12004 IDobler. Keeton. fc Wambsganssl 1200511 . with the 
general conclusion that the amplitude of microlensing 
fluctuations decreases with increasing source size. Thus, 
we naively expect that microlensing fluctuations may 
start out large and decrease in amplitude as the SN ex- 
pands. Another qualitative result is that large sources 
are not very sensitive to microlensing by small stars. 
Therefore we might hope to probe the mass function of 
stars in the lens galaxy by comparing the microlensing 
at different times or source sizes. 

The notion of using multiple source sizes to "map 
out" lens mass distributions on vario us scales has 
been explor ed for quasar l e nsing by iMetcalf et alJ 
i2004D and IDobler fc Keetonl l)2005|) . In particular, 
Dobler & Keetonl plotted curves of magnification versus 
source size and pointed out that multiwavelength ob- 
servations of quasar lenses could measure these curves. 
However, with quasar lenses it is possible to obtain only 
a few discrete points on the curves, say for X-ray, optical, 
and radio wavelengths, and there is no clear mapping be- 
tween the observable quantity (wavelength) and source 
size. With SNe, by contrast, there is a simple mapping 
between the observable quantity (time) and the source 
size via the photosphere expansion rate. Furthermore, 
with good time sampling it would be possible to mea- 
sure the entire magnification versus source size curve. 

Of course, seeing microlensing in SN light curves is con- 
tingent upon having a SN expansion rate that is neither 
too large (which would "wash out" microlensing effects 

3 The sit uation is similar to nano lensing of gamma ray bursts as 
explored by I Walker fc Lewis! <2003l) . 



at early times) nor too small (which would mean that 
the SN remains effectively a point source for its observ- 
able lifetime). In §2 we examine the typical SN expan- 
sion rate in units that are meaningful for microlensing, 
and estimate the optical depth for microlensing by stars 
in both the lens and host galaxies. In §3 we describe 
our microlensing simulations, and in §4 we show sam- 
ple microlensed SNe light curves. In §5 we interpret the 
light curves and consider whether microlensing fluctu- 
tations could be used to extract information about the 
stellar population: the relative abundances of stars and 
smooth dark matter at the image positions, and the mass 
function of stars. We also consider whether microlensing 
will affect time delay measurements in lensed SNe, which 
have b een proposed as a tool for precise measurements 
of H l|Qguri fc Kawanoll2003T) . Finally, in 56 we sum- 
marize our findings and discuss possibilities for future 
theoretical and observational advances. 

2. THEORETICAL BACKGROUND 

2.1. SN model 

For our microlensing calculations, we model the 
source SN as an expanding disk of uniform bright- 
ness, taking as motivation the success of the Expand- 
ing Pho_to^rjbej^_^fethod_(EPM) originally developed 
by iKirshner fc Kwanl l(T974|) . EPM assumes that the 
light emanating from a SN II explosion is due to a 
spherically symmetric, expanding hydrogen shell. The 
method (or closely related variants) has been used by 
sev eral authors to estimate distances to observed SNe 
II iKirshner fc Kwanl 11974 iSchurmann. Arnett. fc FalU 
119791 Branch et alJl!983t lBrancblll987[) . The efficacy of 
EPM rests on two assumptions: the photosphere is spher- 
ical and emits light with a Planck spectrum. Asphericity 
has been measured to be < 30% via spectropolarime- 
try l)Leonard et alll200l|) . and departures from Plankian 
emissi on due to flux dilution has been studied numer- 
ically l|Hershkowitz et all 119861: iHershkowitz fc Wagonerl 
Il987|) . Lastly, for this pilot study we also assume that 
the effects of limb darkening may be neglected. 

2.2. Microlensing scale 

It is important to ask how the size of the SN photo- 
sphere compares to scales relevant for microlensing dur- 
ing the SN's observable lifetime. Given a population of 
stars in the lens galaxy, the relevant length scale for mi- 
crolensing is the Einstein radius 
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associated with the mean stellar mass fh. The distance 
factor D is a function of lens and source redshifts (zi and 
z s respectively) and cosmology. Throughout this paper 
we use zi = 0.32 and z s = 0.53, which represent the mean 
lens and source redshifts of the SLACS sample, and we 
take n M = 0.3, fl A = 0.7, and H = 70 km/s/Mpc. 
With these parameters, equation (JJJ gives 



1/2 



(2) 



To compare this with the size of the SN, we must 
convert the photosphere's expansion rate into units of 
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the Einstein radius. If the edge of the photosphere 
is traveling at a velocity v = 15,00 km/s = 0.05 
c ijEastman. Schmidt. &: Kirshnerlll996j) . this expansion 
rate is 

^ = Zr = 0.0038 R day" 1 , (3) 
dt R 

implying that the size of the SN is less than or compa- 
rable to R for its observable lifetime (t ~ 300 days). In 
other words, the SN has a size that should make it quite 
sensitive to microlensing. 

2.3. Estimates of self-lensing 

Before examining microlensing by stars in an interven- 
ing lens galaxy, it is interesting to consider whether SNe 
(both lensed and unlensed) could be microlensed by stars 
in their own host galaxies. While this effect has yet to be 
identified in known SNe, the large SN catalogs that will 
be produced by SNAP may reveal effects that are rare. 
If we imagine that a SN goes off at the "back" (relative 
to the observer) of its host galaxy, we can estimate the 
self-lensing optical depth 4 by integrating along the line 
of sight through the host galaxy, 
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where 2L is the diameter of the host galaxy, n is the 
number density of stars, and R is again the Einstein ra- 
dius of an individual star. To make a simple toy model, 
we assume that n is constant through the galaxy, that all 
of the stars have the same mass m, and that the galaxy 
can be approximated as a sphere of radius L. 

From the definition of the Einstein radius, equation 10} 
can be written as 



2L 



4nGm didi s 



dl. 



(5) 



where di is the distance to the lens, d s is the distance 
to the source, and di s is separation between the lens and 
the source. Since the distance to the host galaxy will be 
much larger than the galaxy itself, d s « di and di s = £. 
Integrating then yields 
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where we have used the fact that nm = p is the mass 
density of the galaxy, and we have assumed a uniform 
spherical galaxy so p ~ (3m ga i)/ (AirL 3 ) where m ga i is 
the total stellar mass in the galaxy. Note that the self- 
lensing optical depth is independent of the distance to 
the host galaxy. Using typical values for a Milky Way 
type galaxy, m ga i ~ 10 Mq and L ~ 15 kpc, we find 



k* ~ 2 x 10 



(7) 

We conclude that the optical depth is too small, even for 
a SN at the back of the galaxy, to produce a significant 
number of self-lensing events even in the large catalogs 
that will be produced by SNAP. 



4 Following the standard convention t\Qetaj_&^X}striker| I19& 
| Pacvnski|[T98aiKavser. Refsdal. fc Stabelllll98aiWvithe fc Turne 
1200 U) . we define the microlensing optical depth to be the fraction 
of the image plane area covered by the Einstein disks of stars. This 
is equal to the convergence K t , or the surface mass density of stars 
in units of the critical density for lensing. 
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Fig. 1. — Maximum /» as a function of R e i n /R c f{ assuming a 
de Vauco uleurs profile for the lig ht and an isothermal profile for 
the mass. Koopmans ct al. (2006) measure (R a i n /R c s) = 0.52 for 
the SLACS lenses. The maximum surface density in stars at that 
radius is /„ = 0.57. 



2.4. Stellar mass fraction 

The microlensing optical depth is characterized 
by the convergence k*, or the surface mass den- 
sity of star s in units of the cri ti cal densi t y for 
lensing (see iVietri fc Ostrikerl 11983 iPacvnskil Il986t 



Kavser. Refsdal. fc Stabell IT98rT IWvithe & Tumerl 
200lfT It is convenient to define the fraction of the 
surface density in stars, /* = K*/Ktot, at the image 
positions. To make a rough estimate of /*, we model 
the lens galaxy's total mass distribution as a Singular 
Isothermal Sphere (see IKoopmans et al.H2006|) . 
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where i? e ; n is the Einstein radius. We model the stellar 
mass distribution with a dc Vaucouleurs model, 
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where i? c ff is the effective radius, and k — 7.67. If we 
set the stellar mass-to-light ratio equal to its maximum 
possible value such that KDeV ^sis f° r all r, we can 
find the maximum value of the normalization constant A. 
With this maximum value of A we can obtain an upper 
bound on /* = KDeV / ksis as a function of radius. 

Figure ^ shows the upper limit on /» at the Einstein 
radius, as a function of R e - m /R e ff. The average Ein- 
stein radius of the SLA CS lenses is {Rem/ Reft) = 0.52 
l)Koonmans et al. ]|2p06^ , and this corresponds to a max- 
imum stellar mass fraction of /* = 0.57. Also inter- 
esting to consider is the integrated stellar mass fraction 
within the Einstein radius. We find the maximum value 
of this integrated stellar mass fractio n to be 0.83 for 
(R^iaJRcs) — 0.52. For comparison, IKoopmans et alJ 
( 2006) find the mean integrated stellar mass fraction to 
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be 0.75 for the SLACS lenses, suggesting that that the 
actuall stellar mass density is ^90% of its maximum pos- 
sible value. This would imply that /» = 0.46 at the Ein- 
stein radius (again for (R e i a /Res) = 0.52). Therefore, we 
consider /, « 0.5 to be a reasonable estimate of the stel- 
lar mass fraction near the image locations of the SLACS 
lenses 

3. MICROLENSING SIMULATIONS 

We model the population of stars using a Salpeter 
mass function, dn/dra oc m~ 2 35 , spanning the mass 
range from m\ to m,2- We find it convenient to spec- 
ify the mass range by setting the mean mass to be 
fh = Mq and quantifying the dynamic range with q = 
Tn\jm%. Several authors have shown that microlcnsing 
magnific ation distributions are insensitive to the mean 
mass fh ilWvithe fe Turner! 1200 H iSchechter et~ai] 120041 
iMortonson et alJ l2004fl . For microlensing of SN light 
curves, the mean mass does set an overall time scale: 
as fh increases, the time it takes for the source to reach 
a given fraction of R scales as m 1 / 2 . 

The properties of the microlensing fluctuations should 
be dependent on the stellar population parameters 
(namely /» and q). They also depend on the particular 
configuration of stars, because as shown in Figure [5] dif- 
ferent realizations of stellar populations that are statisti- 
cally equivalent yield very different magnification maps 
in the source plane. Note that these magnification maps 
are 2.5.R on a side, so a SN will expand to cover roughly 
half of a map during its observable lifetime. 

We compute microlensing magnifications using a new 
ray-shooting code. The code uses a simplified version 
of a tree algorithm to efficiently compute deflection an- 
gles for a populati on of stars. Where as the ray-shooting 
algorithm used bv lWambseanssI l)1999f) is designed to pro- 
vide magnifications for many source positions (yielding 
magnification maps), our new algorithm is optimized for 
handling a wide range of source sizes. Details of the code 
will be presented elsewhere (Dobler & Keeton, in prep.). 

4. RESULTS 

In our first set of simulations we set k = 7 = 0.45, 
where the total surface mass density, k s is 

the surface density of a smooth dark matter component, 
and 7 is the external shear. We choose k = 7 = 0.45 
for illustrative purposes as these are typical values for 
a positive parity image near the Einstein radius of an 
isothermal lens galaxy. In the absence of microlensing 
(i.e., if there were only "macrolensing" due to the lens 
galaxy), this would yield a magnification of /io = 10. 

Initially we isolate microlensing fluctuations by plot- 
ting the difference light curve, 



AM (t) = -2.5 log 
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where S(t) and I(t) are the intrinsic and observed bright- 
nesses of the SN. This definition is chosen so that AM = 
in the absence of microlensing. In general we expect 
AM(i) — ^ as t — > 00 since large sources are insensi- 
tive to microlensing. Figure |3| shows a sample difference 
ligth curve for one stellar realization with a uniform mass 
function (q — 1) and 100% of the surface density in stars 
(/* = 1). Over the observable lifetime of the SN, this 
stellar configuration creates an offset of up to 0.75 mag, 



and fluctuations at the level of 0.1-0.2 mag. These effects 
are well within the realm of detectability. 

We now combine difference light curves with data for 
the observed SN 2004ej from the Carnegie Supernova 
Project (CSP) 5 . SN 2004ej is a classic example of a SN 
II-P (plateau) in which the light curve falls off slowly 
at early times, drops suddenly by almost 2 mag, and 
then decreases slowly thereafter. Figure 0] shows four ex- 
amples of microlensing effects. Depending on the spe- 
cific stellar configuration, the SN light curve can be 
changed either subtly or dramatically. In particular, 
panel (a) shows a large offset due to microlensing, plus 
small-amplitude fluctuations. Panel (b) shows not only 
an offset of ~1 mag, but also a pronounced peak that 
is not present in the intrinsic light curve. The ability 
of microlensing to induce such large-scale features may 
limit the precision with which one can "match up" light 
curves from different lensed images to measure precise 
lensing time delays (see H5.2|l . Panel (c) shows that the 
microlensing fluctuations can be more subtle with only 
small fluctuations about the intrinsic light curve, while 
panel (d) shows that microlensing can substantially alter 
the decay rate of the SN light curve. 

The lensed image configurations for various epochs of 
panel (b) arc shown in Figure [5j Even when the SN 
is small there are multiple micro-images, and no sin- 
gle image is dominant. As the SN expands, more and 
more small micro-images appear, and a few larger micro- 
images near the center begin to merge. The distribu- 
tion of micro-images is elongated along the shear axis 
(the horizontal axis in Fig. |5J). Even at t — 200 days 
(i? src — 0.76i? c in), by which time several merged im- 
ages have formed, there is still no clear "primary" micro- 
image. The complexity of the image configurations is 
striking and attests the richness of the phenomenon. 

Figure El shows difference light curves for various ran- 
dom star populations (all with q = /„ = 1). It is clear 
that all of the realizations produce offsets, which can be 
as large as 1-1.5 mag. The amplitude of the smaller scale 
fluctuations depends on the timescale. That is, there are 
occasionally very rapid variations at the level of ~0.2 
mag, while longer time scale variations can be as large as 
~0. 5-1.0 mag. To test if these trends are specific to these 
particular stellar population parameters, Figure shows 
the effects of varying q and We show 10 realizations 
of the stellar population for each parameter combination. 
It is clear that microlensing fluctuations are very generic. 
While at first glance it is not obvious whether the realiza- 
tions differ substantially from one panel to another, each 
individual curve has an enormous amount of structure. 
It would seem that a lot of information about the stellar 
population is hidden in the light curves; the challenge is 
figuring out how to extract it (see M5.ll) . 

Figure0also shows that there is not a single realization 
that does not produce fluctuations over the lifetime of the 
SN. To quantify this point, we define 

A max = max(AM(f)) - min(AAf(i)) (11) 

to measure the overall change in any difference curve. 
Figure |H1 shows the probability distribution for A 



max ■ 



5 http://cspl.lco.cl/-cspuserl/CSP.html 

6 We have combined all stellar population parameters and con- 
figurations. While this does wash out detailed information that 
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Fig. 2. — Sample magnificati on maps for four dif ferent realizations of stellar populations with f* = q = 1.0 and m = Mq. (These maps 
were generated with the code by Wambseanss ( 1999), although we actually use our own microlensing code to compute SN light curves.) The 
side length of each map is 2.5R, so an expanding SN will grow to cover approximately one half of a map. The large qualitative differences 
in the magnification maps imply that microlensing fluctuations will be very different from one SN light curve to the next, even for fixed 
stellar population parameters. 



There are two crucial points. First, the A max distribu- 
tion is not peaked at zero, so microlensing fluctuations 
are more likely than not to be observed in SN light curves. 
Second, the total variation in a given SN light curve due 
to microlensing can be quite large. The A max distribu- 
tion is peaked at ~ 0.6 mag with a long tail that extends 
to > 2 mag. Strikingly, the probability that A max is less 
than ~0.1 mag is quite low. From the cumulative distri- 
bution, we estimate that 70% of lensed SN light curves 
will have microlensing fluctuations greater than 0.5 mag, 
and 25% will have fluctuations greater than f mag. We 
have checked that similar conclusions hold for other typ- 
ical values of k and 7 (specifically, those used in *15.2|) . 

may distinguish one set of stellar population parameters from an- 
other, for the moment we only seek to understand how likely it 
is that microlensing fluctuations will be seen in a lensed SN. We 
have examined A m ax distributions for different stellar populations, 
and they are all qualitatitively similar to the distribution shown in 
Figure [S] 



including both positive and negative parity situations. 

The implication is that microlensing of lensed SNe is 
almost a certainty given reasonable stellar populations. 

5. ANALYSIS 

5.1. Extracting parameters of the stellar population 

Comparing the various panels of Figure [7J the eye 
picks out some qualitative features. For example, in the 
q = 0.01 cases there seem to be more regions where the 
difference light curves are fairly flat than in the q = 1.0 
cases. To quantify this feature, we compute the time- 
weighted derivative (TWD) of the difference light curves, 

d 

TWD = -t X —AM, (12) 
dt K ' 

where the negative sign ensures that TWD is positive for 
increasing brightness (i.e., decreasing magnitude). It is 
useful to explain the motivation further. When the mass 
range is broad (q is small), there are a few high-mass 
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Fig. 3. — A sample difference light curve for SN microlensing 
with a uniform mass function (q = 1) and 100% of the surface 
density in stars (/* = 1). Microlensing produces offsets of nearly 1 
mag, and shorter time scale fluctuations as large as 0.1-0.2 mag. 



stars and they tend to be well separated. When the pho- 
tosphere encounters one of these stars there can be an 
abrupt change in the lensing magnification which leads 
to a "step" in the light curve. Conversely, between these 
abrupt events the magnification fluctuations are small. 
A similar argument holds when the stars are diluted by 
a smooth mass component (/* decreases). While differ- 
ent stellar configurations certainly produce different light 
curves, in general the step-type fluctuations are more 
characteristic of low q or low /». The time- weighting 
of the derivative accounts for finite source effects: the 
change in brightness will happen over a longer time scale 
when a larger source comes into contact with a stellar 
caustic. 

The TWD curves are shown in FigurcEl There is strik- 
ing contrast among the different /* and q combinations. 
In general, the q = 0.01 cases do seem more tightly bound 
about TWD = (i.e., flat difference light curves); this 
effect is particularly notable for /* = 0.2. There is an 
even more pronounced effect when comparing different 
/* values at fixed q. Diluting the stellar component with 
smooth dark matter clearly increases the scatter about 
TWD = 0. 

Notice that the TWD curves show large positive spikes, 
but no corresponding negative spikes. This is because 
there is a simple lower bound on the TWD, which can 
be seen as follows. The difference light curve may be 
written as 



AM = -2.5 log 



M07T-R2-C 



(13) 



where Aj is the area of the images. Therefore the TWD 
is 



TWD = 



2.5t 
ln(I0) 



1 dAi 2 dRs 
17 ~di 



R S1C dt 



For linear expansion R s 
plifies and we obtain 

TWD = 



v t so the second term sim- 



2.5 1 dAj 



Ajln(10) dt ln(10) 



(15) 
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In the first term, it is important to understand that 
dAi/dt > 0. This is because as the source expands and 
crosses more and more caustics, images can never disap- 
pear. Distinct images may merge, but they cannot vanish 
because, if a particular source position is covered by the 
photosphere at one time, it is covered at all later times 
as well. 7 This means that the first term is non-negative, 
so we have TWD > —5/ ln(10) = —2.17. In practice, our 
TWD curves do not even dip below TWD ~ — 1, which 
suggests that there may be some other effect bounding 
the TWD. Nevertheless, this argument does explain why 
there are no large negative spikes, only large positive 
spikes. 

To quantify our sense that different stellar populations 
have TWD curves that are more or less tightly bound 
around zero, we compute the dispersion in the TWD dis- 
tribution averaged over all times and all curves for each 
panel in Figure El Figure HTH shows how the TWD dis- 
persion depends on /» for the various values of q. We 
find a clear decrease in the TWD dispersion as /* in- 
creases. There is no obvious trend with q. We conclude 
that the TWD contains key information about the rela- 
tive amounts of stellar versus smooth dark matter mass 
densities. Constraints on the stellar mass function will 
be more difficult to extract from the TWD, but perhaps 
other statistical measures can be found that will reveal 
information about the mass function. 

5.2. Time delay uncertainties 

One goal for lensed SNe is to measure lensing time de- 
lays very precisely. We now consider whether microlens- 
ing may make it difficult to "match up" light curves to 
determine the time delay. To consider a realistic lensing 
scenario, we consider the S LACS lens SPSS J 2300+ 0022 
and use the lens model bv iKoonmans et~aT] p00l . We 
place a SN near the center of the source and solve the 
lens equation to find the local convergence and shear at 
the two SN image positions (k = 7 = 0.39 for the pos- 
itive parity image and k — 7 = 0.52 for the negative 
parity image), and also the time delay between the im- 
ages (Ato = 16.5 days). We then generate mock light 
curves starting from two observed cases: SN 2004ej and 
SN 2004ex from the CSP. We seek to focus on the uncer- 
tainties induced by microlensing, so we limit the uncer- 
tainties due to finite sampling by artificially boosting the 
time sampling to ^4 hours with linear interpolation. We 
generate microlensing difference curves using = 0.5 
(see tl2.4JI . a uniform mass function (q = 1.0), and then 
add them to the intrinsic light curves to obtain different 
realizations of lensed SN light curves. We create 20 ran- 
dom realizations for each image, which can be used to 
construct 400 distinct "observations." 

For each "observation," we then attempt to recover the 
time delay by matching up the light curves from the two 

7 The surface brightness of the photosphere may change with 
time , but tha t does not affect our argument because in equations 
1131 and 1141 we only need to specify the area of the images; the 
source surface brightness has factored out. 
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Fig. 4. — Sample light curves for SN 2004ej, both with (dashed) and without (solid) microlensing fluctuations. The intrinsic light curve 
data are taken from the Carnegie Supernova Project. The four panels represent different stellar configurations which lead to numerous 
qualitative changes in the observed light curve. Microlensing can produce (a) a large offset, (b) non-intrinsic features such as peaks, (c) 
small fluctuations, or (d) changes to the decay rate. 



images. Specifically, we minimize the goodness of fit 
1 N 1 

(16) 

Here D + and D are the positive and negative par- 
ity SN light curves (in magnitudes), respectively. The 
normalization constant k is required because the two 
images have different macro-magnifications from the 
main lens galaxy; in the absence of microlensing, k = 
— 2.51og(/x+//x_). Also, N is the number of data points 
for which the shifted light curves overlap. We normal- 
ize x 2 by -/V in order to avoid a scenario in which the 
optimization routine takes \ 2 to zero simply by shifting 
the light curves so they do not overlap and there are no 
terms in the sum. For simplicitly, we assume that the 
photometric uncertainties Cj are Gaussian and constant 
in time. 

Figure ^2 shows sample observations, together with 
histograms of the recovered time delays. The time delay 



histograms are broad: the dispersions are a — 1.5 days 
for SN 2004ej and 3.6 days for SN 2004ex, implying that 
microlensing will make it difficult to measure time delays 
more precise than a few days. Furthermore, the At his- 
tograms are peaked not at the true time delay but rather 
at a higher value, and are significantly skewed. 

The differences between the results for the two differ- 
ent intrinsic SNe can be understood by comparing the 
intrinsic light curves. SN 2004ej has a very distinct, 
large amplitude intrinsic feature at t w 340 days. The 
light curve decreases rapidly by ^2 mag, which is larger 
than typical microlensing effects (see Fig. [7J). In this 
case, microlensing can perturb the time delay analysis 
but cannot fundamentally disrupt the ability to line up 
the light curves and measure the time delay. Conversely, 
SN 2004ex has intrinsic variability over longer time scales 
and, more importantly, of smaller amplitude. Microlens- 
ing introduces fluctuations with similar time scales and 
amplitudes, which makes it considerably more difficult 
to match up the light curves of the two images. 
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Fig. 5. — Micro-image configurations at various epochs for q = /* = 1.0. As the SN grows in size, micro-images appear and merge. The 
source size of the photosphere for these image configurations is Rs rc = 0.019, 0.038, 0.076, 0.19, 0.38, and 0.76 i? e in f° r times t = 5, 10, 
20, 50, 100, and 200 days respectively. The photometric effects on the intrinsic SN light curve for this particular stellar configuration are 
shown in panel (b) of Figure IH 
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Fig. 6. — Difference light curves for fixed /* = q = 1.0 but 
varying stellar configurations. Short time scale fluctuations ~0.2 
mag are quite common, while longer time scale variations can be 
as large as 0.5—1.0 mag. 



We consider whether it is possible to improve the time 
delay measurements by "masking out" portions of the 



light curve when the source is more susceptible to mi- 
crolensing. Specifically we repeat the time delay analysis 
while ignoring the first 25% of the light curve when the 
source is small and microlensing fluctuations are largest, 
and the last 25% when the source is becoming faint. For 
SN 2004ej this effectively isolates the large intrinsic drop 
around t ~ 340 days, and leads to a smaller time delay 
dispersion of a = 0.6 days. In other words, if the light 
curve has a strong intrinsic feature, it can pay to focus 
on this feature for measuring the time delay. For SN 
2004 ex, by contrast, masking the first and last parts of 
the light curve leads to a larger dispersion of a = 4.40 
days. This is not surprising since, with no dominant in- 
trinsic feature, masking out part of the data only serves 
to decrease the number of data points. 

In other words, the extent to which microlensing cor- 
rupts time delay measurements will depend on whether 
there are strong intrinsic features in the SN light curve. 
Of course, for a real lensed SN the intrinsic light curve 
will not be known, and the challenge will be to determine 
which observed features are intrinsic and which are due 
to microlensing effects. Nevertheless, we can conclude 
in general that microlensing effects must be taken into 
account when attempting to extract high precision time 
delays from multiply imaged SNe. 

6. CONCLUSIONS 

In anticipation of the discovery of lensed SNe in the 
near future with SLACS and SNAP, we have undertaken 
a pilot study of the fluctuations in the observed light 
curves due to microlensing by stars in the lens galaxy. 
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Fig. 7. — Difference light curves for various q and /*. Ten realizations of the stellar configuration arc shown for each case. Significant 
microlensing fluctuations are seen in every realization and for every combination of q and /*, implying that microlensing will be quite 
common in lensed SNe. 



Our main findings can be summarized as follows. 

• The microlensing optical depth and the angular size 
of the source SN virtually ensure that microlens- 
ing fluctuations will be present in lensed SN light 
curves. 

• Microlensing can produce an overall brightness off- 
set of ~1 mag, short time scale brightness fluctua- 



tions of ~0. 1-0.2 mag, and longer time scale fluc- 
tuations of ~0. 5-1.0 mag. 

• We estimate that ^70% of lensed SN light curves 
will have microlensing fluctuations larger than 0.5 
mag, and ~25% will have fluctuations larger than 1 
mag. The probability that microlensing effects will 
not be present in any given light curve is < 5%. 
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Fig. 8. — The probability density function (left) and cumulative distribution function (right) for A max , averaged over all stellar population 
parameters and configurations. The PDF is peaked at 0.6 mag with a long tail extending to > 2 mag. The CDF shows that roughly 70% 
of observed SN light curves will exhibit microlensing fluctuations > 0.5 mag, and 25% will show fluctuations > f mag. 



• Microlensing fluctuations can generate dramatic 
and easily observable changes to the intrinsic SN 
light curve. These features will complicate efforts 
to use time delays to measure the Hubble constant 
to high precision. Without detailed knowledge of 
the intrinsic structure of the light curve, time de- 
lays accurate to better than a few days will be dif- 
ficult to obtain. 

We should say that if microlensing fluctuations are not 
observed in lensed SN light curves, there are two possible 
interpretations: either the ratio of stellar mass density 
to smooth dark matter mass density (at the image po- 
sitions) is so low that the microlensing optical depth is 
near zero; or the SN photosphere undergoes an initial pe- 
riod of very rapid expansion (v 3> 15,000 km/s at early 
times) . 

The microlensing fluctuations depend on the stellar 
mass fraction and mass function in the lens galaxy, al- 
though there are such large difference between different 
random star configurations that it may be challenging to 
extract the stellar population parameters. We have intro- 
duced the time-weighted derivative (TWD) of the light 
curve as a quantity that may be useful for this analy- 
sis. The TWD is relevant because of the correlation be- 
tween fluctuation time scales, stellar Einstein radii, and 
the photosphere expansion velocity. We find that TWD 
statistics may be most useful for determining the rela- 
tive abundances of stars and dark matter at the image 



positions. More generally, the amount of structure in the 
light curves is a sign that more effort should be exerted 
to figure out how to analyze lensed SN light curves to 
determine properties of the stellar population. 

Once a lensed SN is typed via spectroscopy or the late- 
time shape of the photometric light curve, it may be 
possible to model the intrinsic light curve. The efficacy 
of EPM implies that the source geometry at different 
epochs in the observations will be fairly well known, and 
so a light curve fit ting method (such as that described 
by Kochanek 2004) may ultimately be the best analysis. 
The fitting method must be extended to efficiently incor- 
porate expanding sources as opposed to moving sources; 
however the basic idea of finding a family of possible stel- 
lar configurations and using a Bayesian likelihood anal- 
ysis to extract mass function parameters will remain the 
same. 
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HST-AR-10668 from the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities 
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Fig. 10. — Dispersion a in the TWD distributions, plotted 
versus /* for the three values of q. There is a clear decrease of a as 
/* is increased from 0.2 to 1.0. Thus, microlensing fluctuations in 
SN light curves may be useful in extracting the relative amounts 
of stellar and smooth dark matter at the positions of the lensed 
images. 
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Fig. 11. — Illustration of the effects of microlensing on time delay measurements. We create a mock lensed SN using two different 
intrinsic SN light curve s (SN 2004ej and SN 2004ex, both from the CSP) and lens/source parameters from the SLACS lens J2300+0022 
IKoopm ans ct al. 2006). Left column: Histograms of the recovered time delays for 400 mock observations. The true time delay Aip = 16.54 
days is marked with a vertical line. The histograms are broad (dispersion cr = 1.48 days and 3.55 days for SN 2004ej and SN 2004ex 
respectively) and skewed with peaks offset from AtQ. The implication is that microlensing effects will render high precision time delays 
difficult to measure from multiply imaged SNe. Center column: The intrinsic light curves for the two SNe. The leading and trailing light 
curves correspond to the positive and negative parity images, respectively. Right column: The positive and negative parity light curves, 
including microlensing fluctuations, for one sample mock observation. 



